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Equilibrium  potential  hysteresis  of  electrode  materials  refers  to  the  phenomenon  that  discharge  equi¬ 
librium  potential  is  lower  than  charge  equilibrium  potential.  It  is  induced  by  the  strain  accommodation 
energy  of  phase  transformation.  Taking  LiFeP04  as  an  example,  the  equilibrium  potential  hysteresis  and 
accommodation  energy  of  two  LiFeP04  with  different  particle  sizes  were  characterized  using  galvano¬ 
static  intermittent  titration  technique  (GITT),  cyclic  voltammetry  (CV),  and  electrochemical  impedance 
spectroscopy  (EIS)  methods.  The  LiFeP04  sample  with  40  nm  particle  size  has  lower  accommodation 
energy  than  that  of  lOOnm  LiFeP04,  which  is  due  to  a  narrower  miscibility  gap  and  smaller  potential 
hysteresis  of  40 nm  LiFeP04  comparing  to  lOOnm  LiFeP04.  The  phase  transformation  of  LiFeP04  occurs 
only  at  potential  levels  above  delithiation  (charge)  equilibrium  potential  and  below  lithiation  (discharge) 
equilibrium  potential.  No  phase  transformation  occurs  when  charge/discharge  is  limited  in  the  potential 
hysteresis  range,  i.e.  LiFeP04  behaves  like  a  solid  solution  in  the  potential  hysteresis  range.  A  reliable  EIS 
can  only  be  obtained  when  it  is  performed  at  the  center  of  potential  hysteresis  and  the  potential  ampli¬ 
tude  does  not  exceed  the  range  of  equilibrium  potential  hysteresis.  The  analysis  on  strain  accommodation 
and  potential  hysteresis  of  LiFeP04  cathodes  during  Li  ion  insertion/extraction  is  also  suitable  for  other 
phase  transformation  electrodes. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

LiFeP04  has  been  considered  to  be  a  promising  candidate  mate¬ 
rial  for  the  cathode  of  high-power,  safe,  low-costand  long-life  Li-ion 
batteries  required  for  hybrid  electric  vehicles  and  renewable  energy 
storage.  The  most  striking  characteristic  of  this  material  is  a  flat 
charge/discharge  profile  at  the  potential  of  ~3.4  V  versus  Li/Li+  due 
to  the  phase  transformation.  However,  pure  LiFeP04  suffered  from 
poor  rate  capability  due  to  low  intrinsic  electronic/ionic  conduc¬ 
tivity  and  slow  phase  transformation.  Extensive  research  has  been 
carried  out  to  overcome  the  poor  reaction  kinetics  of  LiFeP04.  The 
most  effective  strategies  include  (i)  to  downsize  the  LiFeP04  par¬ 
ticles  to  nano-scale  [1]  and  (ii)  to  coat  electronically  conductive 
carbon  [2]  or  an  ionically  conductive  phase  [3]  onto  the  LiFeP04 
particle  surface.  However,  only  few  studies  have  been  focused  on 
the  phase  transformation  [4-7].  Since  most  capacity  of  LiFeP04 
comes  from  a  first-order  phase  transformation,  an  understanding 
of  the  kinetics  of  the  phase  transformation  is  critical  to  achieving 
a  high-power  density  for  LiFeP04,  especially  for  the  carbon-coated 
nano-LiFeP04  where  the  electronic  and  ionic  transport  have  been 
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greatly  improved  and  the  slow  phase  transformation  may  be  a  lim¬ 
iting  step  for  the  rate  performance. 

Recent  research  demonstrated  the  phase  transformation  in 
LiFeP04  is  strongly  controlled  by  the  strain  accommodation  energy 
induced  from  the  volume  difference  between  the  lithiated  and 
delithiated  phases  [4].  The  misfit  strain  between  the  triphylite 
phase  (Lii_yFeP04)  and  the  heterosite  phase  (Li*FeP04)  in  nano- 
LiFeP04  particles  was  detected  using  X-ray  diffraction  (XRD),  [4] 
and  dislocations  induced  by  misfit  strain  were  observed  in  plate¬ 
like  LiFeP04  particles  using  transmission  electron  microscopy 
(TEM)  [8].  Such  a  strain-induced  energy  penalty  can  decrease  the 
driving  energy  for  phase  transformation,  [9]  resulting  in  decrease 
of  the  phase  transformation  rate  [10].  The  accommodation  energy 
decreases  the  lithiation  equilibrium  potential  [9]  but  increases  the 
delithiation  potential,  resulting  in  an  equilibrium  potential  hystere¬ 
sis  [11].  It  was  reported  that  shrinking  the  LiFeP04  particle  size 
can  narrow  the  miscibility  gap  between  the  lithiated  and  delithi¬ 
ated  phases  and  decrease  the  equilibrium  potential  hysteresis  [4]. 
The  fast  phase  transformation  may  be  one  of  the  reasons  for  the 
excellent  rate  performance  of  carbon-coated  nano-LiFeP04. 

Although  the  equilibrium  potential  hysteresis  is  a  critical  prop¬ 
erty  of  phase  transformation  electrodes,  the  characteristics  of 
the  equilibrium  potential  hysteresis  of  LiFeP04  cathodes  have 
never  been  investigated.  In  this  paper,  the  potential  hysteresis  of 
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two  LiFeP04  samples  with  different  particle  sizes  was  analyzed 
using  galvanostatic  intermittent  titration  technique  (GITT),  cyclic 
voltammetry  (CV)  and  electrochemical  impedance  spectroscopy 
(EIS). 


2.  Experimental  methods 


Two  LiFeP04  samples  were  provided  by  a  commercial  supplier. 
Transmission  electron  microscopy  (TEM)  images  for  both  samples 
were  recorded  on  a  Flitachi  FIF-2000  FEG-TEM,  shown  in  Fig.  1. 
The  LiFeP04  samples  consisted  of  aggregates  of  secondary  nano¬ 
scale  crystal  particles  with  average  sizes  of  40 nm  (Fig.  la)  and 
lOOnm  (Fig.  lb).  The  thin  carbon  coating  (2-4 nm)  and  crystalline 
single  phases  of  LiFeP04  can  be  clearly  observed  in  the  high  res¬ 
olution  TEM  images  (Fig.  lc),  in  which  the  boundaries  between 
LiFeP04  crystal  and  carbon  coating  have  been  marked.  Particle  size 
distributions  of  two  LiFeP04  samples  were  obtained  by  randomly 
measuring  50  particles  from  the  TEM  images.  The  mean  particle  size 
of  two  LiFeP04  samples  was  calculated  by  the  following  equation: 


dm  — 


EM 

Eim 


(i) 


where  n,-  is  the  number  of  particles  with  diameter  d/.  Two  LiFeP04 
samples  showed  very  narrow  particle  size  distribution. 

All  electrochemical  tests  in  this  paper  were  performed 
using  three-electrode  pouch  cells.  The  LiFeP04  cathodes  con¬ 
taining  82wt%  active  material,  10wt%  carbon  black  and  8wt% 
poly(vinylidesfluoride)  (PVDF)  in  l-methyl-2-pyrrolidinone  (NMP) 
solvent  were  prepared  by  the  slurry  coating  method.  The  formula¬ 
tion  was  coated  onto  aluminum  foil  with  a  loading  of  2.5  mg  cm-2  of 
LiFeP04.  The  three-electrode  pouch  cells  were  assembled  in  an  Ar- 
filled  glove  box,  using  lithium  metal  as  the  reference  and  counter 
electrodes,  1  M  LiPF6  in  1 : 1 : 1 :3  ethylene  carbonate  (EC):dimethyl 
carbonate  (DMC):diethyl  carbonate  (DEC):ethylmethyl  carbonate 
(EMC)  as  electrolyte  and  Celgard  3501  as  the  separator. 

To  measure  the  relationship  between  equilibrium  potential 
(fully  relaxed  open-circuit-potential  (OCP))  and  state  of  charge 
(SOC)  or  state  of  discharge  (SOD),  cells  were  tested  using  a  gal¬ 
vanostatic  intermittent  titration  technique  (GITT).  All  the  GITT  tests 
have  been  performed  using  Arbin  test  station.  GITT  measurements 
consisting  of  a  series  of  current  pulses  were  applied  to  the  three- 
electrode  pouch  cells  at  a  low  current  of  less  than  0.022  C  for  2  h, 
each  followed  by  a  1 6  h  rest  process.  The  OCP  at  the  end  of  the  1 6th 
hour  is  considered  to  be  the  equilibrium  potential.  The  potential 
decay  during  the  rest  process  includes  both  relaxation  and  self¬ 
discharge.  LiFeP04  electrodes  cannot  be  fully  relaxed  to  equilibrium 
if  the  relaxation  time  is  too  short.  However,  a  relaxation  time  that 
is  too  long  will  greatly  increase  the  testing  time,  resulting  in  notice¬ 
able  self-discharge.  The  potential  decay  will  be  mainly  controlled 
by  self-discharge  if  relaxation  time  is  much  longer  than  16h,  as 
evidenced  from  the  capacity  difference  between  charge  GITT  and 
discharge  GITT.  The  relaxation  time  of  16  h  was  selected  to  allow 
full  relaxation  of  OCP  and  to  minimize  the  self-discharge  of  LiFeP04 
during  the  test.  The  potentials  of  2.2  and  4.2  V  were  used  as  the  low 
and  high  cutoff  voltages  in  GITT  tests. 

To  characterize  the  equilibrium  potential  hysteresis  of  LiFeP04 
in  the  phase  transformation  region,  cyclic  voltammetry  (CV)  at  a 
scan  rate  of  0.05  mV  s-1  was  performed  on  electrodes  at  the  mid¬ 
dle  of  potential  hysteresis  between  50%  SOC  and  50%  SOD  using 
different  scan  potential  amplitudes.  A  low  scan  rate  of  0.05  mV  s-1 
is  used  here  to  decrease  the  influence  of  lithium  ion  diffusion  on  the 
current  but  detect  the  phase  transformation  current  [12].  Similarly, 
electrochemical  impedance  spectroscopy  (EIS)  was  also  applied  on 
the  LiFeP04  cathodes  at  50%  SOC,  50%  SOD  and  at  the  middle  of 
potential  hysteresis  between  50%  SOC  and  50%  SOD.  Before  CV  and 


EIS  tests,  the  electrodes  were  charged/discharged  between  2.2  and 
4.2  V  at  0.1  C  current  for  two  cycles.  Both  CV  and  EIS  were  tested  on 
Solartron  1260/1287. 

3.  Results  and  discussion 

3.1.  Equilibrium  potential  hysteresis  ofLiFeP04  induced  by 
accommodation  energy 

Delithiation/lithiation  of  LiFeP04  is  usually  accompanied  by 
phase  transformation  between  heterosite  phase  (a  phase)  and  tri- 
phylite  phase  (|3  phase)  [4].  The  formation  of  potential  hysteresis 
during  delithiation  and  lithiation  of  LiFeP04  was  discussed  in  our 
publication  [  1 3  ].  In  brief,  during  Li  ion  insertion,  the  lattice  parame¬ 
ters  of  the  LixFeP04  vary  with  the  lithium  stoichiometry  and  phase 
transformation,  since  the  newly  formed  phase  has  a  different  molar 
volume  with  that  of  the  matrix  [4,14].  To  accommodate  the  vol¬ 
ume  difference,  stress  and  strain  will  be  generated  during  Li  ion 
insertion  into  FeP04.  For  a  purely  elastic  matrix,  elastic  energy 
will  be  stored  in  the  matrix  and  the  newly  formed  phase  during 
phase  transformation  [15].  For  a  plastic  matrix  with  a  low  yield 
strength,  the  elastic  energy  stored  in  the  matrix  will  be  relaxed 
and  spent  in  creating  plastic  deformation  [16].  Normally,  the  elec¬ 
trode  materials  are  elastic-plastic  solid,  the  transformation  energy 
is  accommodated  by  both  elastic  and  plastic  energy  [15].  The  plas¬ 
tic  deformation  in  LiFeP04  during  phase  transformation  has  been 
evidenced  from  the  existence  of  dislocations  and  fracture  [8,17]. 
The  elastic-plastic  accommodation  energy  is  an  energy  barrier 
against  lithiation  and  delithiation  [9].  It  increases  the  free  energy 
of  the  material  and  decreases  the  lithiation  equilibrium  poten¬ 
tial.  The  decrease  of  lithiation  equilibrium  (open-circuit)  potential 
AF lithiation  can  be  determined  from  [15,18]. 

^ relastic  ,  ^ ^plastic 

Ac  n  n  L7 lithiation  '  lithiation 

lithiation  =  ho  ~  tde  =  - ~p -  M) 

where  E0  is  the  theoretical  (accommodation- free)  equilibrium 
potential,  Ede  is  the  discharge  equilibrium  potential  after  consid¬ 
eration  of  the  accommodation  process  during  Li  ion  insertion, 
A(SL  and  AC^L  are  elastic  and  plastic  accommodation 
energy  during  lithiation  process,  n  is  the  number  of  electrons  passed 
per  atom  of  host  material  reacted  and  F  is  the  Faraday’s  constant. 
The  accommodation  energy  increases  with  the  depth  of  discharge 
process,  which  is  similar  to  the  strain-hardening  phenomenon 
observed  in  phase  transformation  of  metal  hydrides  [19]. 

During  Li  ion  extraction,  a  reverse  phase  transformation  will 
occur.  The  volume  changes  during  the  Li  ion  extraction  will  also 
be  accommodated  by  the  elastic-plastic  process  outlined  above. 
Different  from  the  lithiation  process,  the  accommodation  energy 
during  Li  ion  extraction  increases  the  electrode  potential.  The 
potential  increase  during  the  charge  (delithiation)  process  can  be 
calculated  as 

^ relastic  ,  ^rplastic 

\  r  r  r  VJ delithiation  '  VJ delithiation 

delithiation  =  bce  -  t0  =  - — -  (JJ 

where  Ece  is  the  charge  equilibrium  potential  after  consideration 
of  the  accommodation  process,  AC^Mn  and  AGj£Lon  are 
elastic  and  plastic  accommodation  energy  during  delithiation  pro¬ 
cess.  The  accommodation  processes  during  Li  ion  insertion  and 
extraction  dissipate  energy  to  the  surroundings,  resulting  in  an 
open-circuit-potential  hysteresis  (AEUthiation  +  A EdeUthiation).  There¬ 
fore,  the  equilibrium  (open-circuit)  potential  hysteresis  during  Li 
ion  insertion  and  extraction  can  be  calculated  from  the  accommo¬ 
dation  energy.  From  the  accommodation  energy  of  Li-Sn  alloy,  Hirai 
et  al.  [9]  calculated  the  lithiation  equilibrium  potential  of  Sn  using 
Eq.  (2),  which  is  in  excellent  agreement  with  measured  data.  Eqs.  (2) 
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Fig.  1.  TEM  images  for  (a)  LiFeP04  with  secondary  particle  size  ~40nm  and  (b)  LiFeP04  with  secondary  particle  size  ~100nm.  (c)  The  carbon  coating  layer  on  a  crystal 
LiFeP04  particle  is  marked  in  HTEM. 


and  (3)  have  also  been  applied  in  the  metal-hydrogen  system  for 
determining  the  pressure  hysteresis  of  metal  hydrides  from  their 
accommodation  energy  [11]. 

Equilibrium  potential  hysteresis  of  phase  transformation  elec¬ 
trode  during  lithiation  and  delithiation  can  easily  be  determined 
using  GITT,  however,  the  accommodation  energy  is  very  difficult 
to  obtain  from  materials  properties  especially  for  new  electrode 
materials.  Based  on  the  relationship  between  potential  hysteresis 
and  accommodation  energy  in  Eqs.  (2)  and  (3),  the  accommodation 
energy  can  be  calculated  from  equilibrium  potential  hysteresis. 

The  following  reaction  occurs  during  phase  transformation  from 
heterosite  (Li*FeP04)  to  triphylite  (Li1_yFeP04). 

lithiation 

Li*FeP04  +  (l -y-x)Li+  +  (l  -y-x)e~  =  Li,_yFeP04  (4) 

delithiation 

When  the  accommodation  energy  =  A Gg'“*on  + 

AC jj'thjation )  change  during  lithiation  is  defined  as  joules  per  molar 
o/FeP04,  the  accommodation  energy  for  the  lithiation  process  can 
be  expressed  as 

A G~n  =  O  -  ■ * -yWithiation  (5) 

Since  the  values  of  x,  1  -y  and  AEiithiation  can  be  obtained  from 
the  equilibrium  potential-composition  isotherm  curves,  [15]  the 
accommodation  energy  can  be  calculated  using  Eq.  (4).  Similarly, 
the  accommodation  energy  during  delithiation  can  also  be  obtained 
using 

A  GSSSton  =0-*-  y)F^delmation  (6) 

The  equilibrium  potential  of  LiFeP04  with  secondary  particle 
sizes  of  40  and  lOOnm  were  tested  at  different  levels  of  lithia¬ 
tion  and  delithiation  (Fig.  2).  The  equilibrium  potentials  of  both 
LiFeP04  electrodes  gradually  decreased  with  Li  ion  insertion  but 
increased  with  Li  ion  extraction.  The  variation  in  equilibrium  poten¬ 
tial  of  40  nm  LiFeP04  during  Li  ion  insertion  and  extraction  is  less 
than  that  of  lOOnm  LiFeP04.  Both  LiFeP04  particles  have  a  sim¬ 
ilar  delithiation  equilibrium  potential  within  the  miscibility  gap, 


however,  the  lithiation  equilibrium  potential  of  40  nm  LiFeP04  is 
5-6 mV  higher  than  that  of  lOOnm  LiFeP04,  resulting  in  a  smaller 
potential  hysteresis  (8-10  mV)  for  40  nm  LiFeP04,  compared  to 
13-16 mV  for  lOOnm  LiFeP04.  In  addition,  nano-scale  (40 nm) 
LiFeP04  has  a  reduced  miscibility  gap  compared  to  coarser-grained 
(lOOnm)  materials.  The  size-dependent  potential  hysteresis  and 
miscibility  gap  shown  in  Fig.  2  are  well  in  agreement  with  the 
observations  from  other  researchers  [4,5,20]. 

The  accommodation  energies  of  LiFeP04  with  different  particle 
sizes  during  the  first  cycle  of  lithiation  and  delithiation  in  the  phase 
transformation  region  were  calculated  from  equilibrium  potential 
hysteresis  using  Eqs.  (5)  and  (6)  and  shown  in  Fig.  3.  The  theo¬ 
retical  equilibrium  potential  E0  is  defined  as  the  average  value  of 
potentials,  at  which  the  charge  and  discharge  phase  transformation 
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Fig.  3.  Discharge  and  charge  accommodation  energy  for  LiFeP04  with  different  par¬ 
ticle  sizes. 

initiated.  Compared  to  100  nm  LiFeP04, 40  nm  LiFeP04  has  a  lower 
accommodation  energy  due  to  a  narrow  potential  hysteresis  and 
a  reduced  miscibility  gap  in  nano-size  LiFeP04  [5].  Both  the  short 
diffusion  length  and  a  low  accommodation  energy  in  nano-LiFeP04 
accelerate  the  phase  transformation  rate,  resulting  in  high  rate  per¬ 
formance  [4].  The  accommodation  energy  of  both  LiFeP04  samples 
increased  with  SOD  and  SOC.  At  Li0.5FePO4,  the  discharge  (lithia- 
tion)  accommodation  energy  is  larger  than  the  charge  (delithiation) 
accommodation  energy.  The  reduced  accommodation  energy  in 
charge  compared  to  discharge  may  be  one  of  the  reasons  for  higher 
charge  rate  capability  of  LiFeP04  than  that  for  discharge  [21  ]. 

It  has  been  reported  that  dislocations  and  other  defects  were 
generated  during  the  first  cycle  of  Li  ion  insertion/extraction  in 
LiFeP04  [8,17].  These  defects  change  the  accommodation  energy 
in  the  following  cycles.  Therefore,  the  accommodation  energy 
in  the  second  charge/discharge  cycle  may  be  different  from  the 
values  in  the  first  cycle.  Fig.  4a  compares  the  equilibrium  poten¬ 
tials  of  lOOnm  LiFeP04  in  the  first  and  second  charge/discharge 
cycles.  The  corresponding  accommodation  energies  in  the  first 
and  second  charge  and  discharge  cycles  are  illustrated  in  Fig.  4b. 
As  shown  in  Fig.  4a,  the  discharge  equilibrium  potential  in  the 
second  charge/discharge  cycle  is  almost  the  same  as  that  in  the 
first  cycle,  but  the  charge  equilibrium  potential  in  the  second 
charge/discharge  cycle  is  lower  than  that  in  the  first  cycle  at  later 
SOC.  The  lower  charge  equilibrium  potential  induced  lower  charge 
accommodation  energy  in  the  second  charge  cycle.  These  behav¬ 
iors  are  most  likely  due  to  (i)  the  introduction  of  dislocations 
and/or  fractures  in  micro-LiFeP04  in  the  first  charge/discharge 
cycle,  which  lowered  the  energy  required  for  phase  growth  for 
the  second  charge  cycle,  and  (ii)  the  volume  expansion/shrinkage 
of  LiFeP04  during  the  first  charge/discharge  cycle,  which  may 
locally  rearrange  the  position  of  LiFeP04  particles  to  accommodate 
the  stress.  The  reduced  accommodation  energy  during  the  initial 
charge/discharge  cycles  is  consistent  with  the  activation  process 
of  LiFeP04  electrodes.  Flowever,  why  these  defects  only  affect  the 
charge  accommodation  energy  but  not  the  discharge  accommoda¬ 
tion  energy  is  not  clear. 

As  shown  in  Fig.  4b,  the  discharge  accommodation  energy  of 
lOOnm  LiFeP04  increased  up  to  750 J mol-1  when  LixFeP04  (a 
phase)  gradually  transformed  into  Li!_yFeP04  (p  phase).  The  maxi¬ 
mum  accommodation  energy  (750 J  mol-1 )  of  LiFeP04  during  phase 
transformation  from  LixFeP04  to  Lii_yFeP04  is  less  than  the  accom¬ 
modation  energy  (972  J  mol-1 )  of  Nb  during  phase  transformation 


Fig.  4.  (a)  GITT  equilibrium  potential  at  different  cycles  and  (b)  charge  and  discharge 
accommodation  energy  at  different  cycles. 

from  Nb  to  NbFI  [15]  and  is  much  smaller  than  the  accommodation 
energy  (9630 J  mol-1 )  of  Sn  during  phase  transformation  from  Sn  to 
Li0.4Sn  [9].  Therefore,  the  exceptional  rate  performance  of  LiFeP04 
may  be  partially  attributed  to  the  lower  accommodation  energy 
during  phase  transformation  compared  to  Sn  anode  for  Li-ion  bat¬ 
tery  and  Nb-hydride  anodes  for  Ni/MH  batteries. 

3.2.  Characterization  of  potential  hysteresis  with  GITT 

To  characterize  the  property  of  equilibrium  potential  hysteresis 
of  LiFeP04,  the  delithiation  equilibrium  potential  of  1 00  nm  LiFeP04 
was  measured  after  being  equilibriately  discharged  to  different  lev¬ 
els  (SOD:  40%,  60%,  80%  and  100%)  using  GITT  and  shown  in  Fig.  5. 
At  a  given  Li  ion  concentration,  the  charge  equilibrium  potential 
increased  with  the  depth  of  the  previous  SOD.  For  example,  the 
electrode  with  composition  Li0.3FePO4>  which  was  previously  com¬ 
pletely  discharged,  has  the  highest  charge  equilibrium  potential, 
while  the  Li0.3FePO4  electrode  that  was  previously  discharged  to 
Li0.4FePO4  has  the  lowest  charge  equilibrium  potential.  Although 
the  equilibrium  potential  of  the  following  charge  process  increased 
with  the  depth  of  previous  SOD,  the  initial  potential  for  the  phase 
transformation  from  p  to  a  phase  during  charge  process  occurs 
at  the  same  value  (see  A-B  line  in  Fig.  5)  no  matter  how  much 
Li  ion  was  inserted  into  FeP04  in  the  previous  discharge.  If  the 
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Fig.  5.  Charge  equilibrium  potential  of  100  nm  LiFeP04  following  SOD  40%,  60%,  80, 
and  100%.  Discharge  and  charge  current:  0.022  C,  discharge  time:  2h,  charge  time: 
0.33  h  for  40%  and  0.67  h  for  60,  80  and  100%.  Rest  time  for  charge  and  discharge: 
8  h. 

charge  equilibrium  potential  lines  previously  discharged  to  SOD 
40%,  60%  and  80%  are  moved  along  the  A  ->  B  direction  to  the  charge 
potential  line  previously  discharged  to  SOD  100%,  four  charge  lines 
will  be  completely  overlapped.  Therefore,  the  charge  equilibrium 
potentials  following  different  SOD  levels  are  solely  attributed  to 
the  fraction  of  newly  formed  a  phase  rather  than  the  fraction  of  a 
phase  left  during  previous  SOD. 

To  illustrate  the  kinetics  of  charge  and  discharge  of  LiFeP04, 
the  potential  responses  during  each  current  pulse  in  the  GITT 
measurement  were  recorded  and  shown  in  Fig.  6.  Due  to  the  res¬ 
olution  limitation  of  Arbin  equipment,  the  potential  jump  due  to 
solid-electrolyte-interphase  resistance  and  charge  transfer  reac¬ 
tion  cannot  be  monitored.  The  polarization  curves  recorded  in 
Fig.  6  only  reflect  potential  change  due  to  Li-ion  diffusion  and 
the  phase  transformation.  As  demonstrated  in  Fig.  6,  both  charge 
and  discharge  overpotentials  increased  with  the  fraction  of  new 
phase.  To  investigate  the  kinetics  of  discharge  at  different  SOD, 
the  GITT  potential-time  profiles  at  the  2nd,  5th  and  10th  discharge 
pulse  were  compared  by  shifting  their  initial  open-circuit  poten¬ 
tials  together,  as  shown  in  Fig.  7.  The  changes  in  equilibrium  OCP 
( AE  in  Fig.  7)  before  and  after  current  pulses  are  almost  the  same 
in  three  current  pulses.  However,  overpotential  (the  potential  dif¬ 
ference  between  fully  relaxed  OCP  and  the  potential  at  the  end  of 
the  current  pulse)  in  Fig.  7  increased  with  SOD.  By  carefully  com¬ 


Fig.  7.  Discharge  and  following  rest  potential  during  GITT  at  the  2nd,  5th  and  10th 
current  pulses. 

paring  the  potential  response  in  each  current  pulse  (Fig.  7),  we 
found  that  the  initial  potential  exponentially  drops,  followed  by  an 
almost  linear  decrease  with  time.  The  initial  potential  drop  is  due 
to  the  Li-ion  diffusion,  and  the  later  slow  decrease  in  potential  is 
attributed  to  the  phase  transformation.  The  diffusion  polarization 
quickly  increases  with  SOD,  but  the  phase  transformation  polar¬ 
ization  only  slightly  changes  with  SOD.  Fig.  7  demonstrates  that 
the  increase  in  overpotential  with  SOD  is  mainly  attributed  to  the 
increase  in  Li-ion  diffusion  length  due  to  the  growth  of  the  (3  phase. 
The  increase  in  diffusion  overpotential  with  phase  transformation 
was  also  observed  in  the  lithiation/delithiation  of  graphite  [22]. 

The  charge  GITT  following  different  SOD  was  used  to  char¬ 
acterize  potential  hysteresis.  As  shown  in  Fig.  6  the  polarization 
curves  are  SOD-independent,  i.e.  the  potential  responses  in  charge 
GITT  following  different  SOD  were  similar  at  the  same  potential 
level  except  for  40%  SOD  in  which  the  charge  time  of  LiFeP04  was 
shortened  from  0.66  h  (for  all  other  SOD)  to  0.33  h  to  get  more  equi¬ 
librium  potential  data.  To  clearly  show  the  potential  response,  the 
discharge  and  charge  GITT  curves  at  60%  SOD  were  enlarged  and 
shown  in  Fig.  8.  The  potential  response  in  each  charge  GITT  curve 
shows  only  exponential  increase  (diffusion  polarization)  without 
linear  polarization  (phase  transformation  polarization)  until  the 
7th  charge  GITT  curve,  where  the  polarization  potential  exceeds 
the  charge  equilibrium  potential  and  phase  transformation  occurs. 
Therefore,  no  phase  transformation  occurs  inside  potential  hys- 


x  in  LixFeP04 


Fig.  6.  Charge  GITT  curves  following  different  SOD. 


Fig.  8.  Charge  and  discharge  GITT  at  60%  SOD. 
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Fig.  9.  CV  performed  at  the  middle  of  voltage  hysteresis  with  different  amplitudes. 


teresis  region.  Although  in  the  potential  hysteresis  region,  LiFeP04 
still  consists  of  two  phases,  it  behaves  like  a  solid  solution.  Here, 
we  need  to  point  out  that  during  charge/discharge,  in  the  poten¬ 
tial  hysteresis  region,  Li  ions  can  penetrate  across  the  interface, 
and  the  interface  between  two  phases  in  LiFeP04  does  not  move, 
which  is  different  from  the  behavior  of  phase  transformation  and 
similar  to  solid  solution.  The  solid  solution  behavior  in  the  poten¬ 
tial  hysteresis  region  is  also  supported  by  the  fact  that  the  charge 
polarization  curves  after  different  SOD  are  similar  to  the  charge 
polarization  curves  after  fully  discharged  (100%  SOD),  shown  in 
Fig.  6,  and  the  fully  discharged  ( 1 00%  SOD)  LiFeP04  is  a  solid  solution 
before  the  charge  potential  reaches  the  charge  equilibrium  poten¬ 
tial.  In  addition,  the  first  charge  overpotential  is  much  smaller  than 
the  previous  discharge  overpotentials.  For  example,  the  overpo¬ 
tential  at  the  10th  discharge  GITT  (lOd  in  Fig.  8)  is  around  6.0  mV 
at  the  discharge  current  of  0.022  C.  However,  the  overpotential  in 
the  first  charge  GITT  (lc  in  Fig.  8)  is  only  3  mV  at  the  same  cur¬ 
rent,  while  the  overpotential  in  the  following  discharge  GITT  (lid 
in  Fig.  8)  is  7  mV.  The  difference  between  the  11th  discharge  and 
the  1st  charge  is  that  phase  transformation  occurs  in  the  11th  dis¬ 
charge,  while  no  phase  transformation  takes  place  in  the  1  st  charge. 
The  phase  transformation  overpotential  at  the  11th  discharge  GITT 
is  around  4.0  mV,  which  is  57%  of  the  total  (phase  transformation 
and  diffusion)  overpotential. 

3.3.  Characterization  of  phase  transformation  potential  hysteresis 
with  CV  and  EIS 

Cycle  voltammetry  (CV)  is  a  common  and  useful  tool  to  dis¬ 
tinguish  the  phase  transformation  from  a  solid  solution  because 
the  current  will  be  enhanced  in  the  phase  transformation  region. 
Therefore,  the  CV  is  used  to  characterize  the  phase  transformation 
behavior  of  LiFeP04  across  the  potential  hysteresis  where  the  phase 
transformation  of  LiFeP04  occurs  only  when  the  applied  potential 
is  lower  than  the  lithiation  equilibrium  potential  or  higher  than 
delithiation  equilibrium  potential.  LiFeP04  in  the  potential  hys¬ 
teresis  range  can  be  considered  as  a  solid  solution  material.  Fig.  9 
shows  the  CV  scans  of  lOOnm  LiFeP04  at  the  center  of  poten¬ 
tial  hysteresis  (marked  in  Fig.  4a)  using  different  scan  potential 
amplitudes  but  a  fixed  scan  rate  of  0.05  mV  s-1.  Before  each  CV 
measurement,  the  lOOnm  LiFeP04  electrode  was  fully  charged  to 
4.2  V  and  then  discharged  to  SOD  50%  followed  by  charging  the 
electrode  to  3.427  V  and  holding  the  electrode  at  this  potential 
for  2  h.  3.427  V  is  the  potential  halfway  between  lithiation  poten¬ 
tial  (3.422  V)  and  delithiation  potential  (3.432  V).  By  changing  the 
amplitude  of  CV  (±2,  ±4,  ±6,  ±8,  ±10,  ±12,  ±15,  ±20,  ±25,  ±30, 
±35  mV),  the  reversibility  of  the  reactions  inside  and  outside  of  the 


Fig.  10.  EIS  plots  measured  at  (a)  OCP  of  50%  SOD  and  50%  SOC,  and  (b)  the  middle 
potential  of  the  potential  hysteresis  with  different  voltage  amplitudes. 


potential  hysteresis  range  was  characterized.  When  the  amplitudes 
of  CV  are  similar  or  less  than  the  half  value  of  potential  hystere¬ 
sis  (6-8  mV),  the  potential  almost  linearly  changed  with  current. 
However,  for  CV  with  a  large  amplitude  (>12  mV),  a  rapid  increase 
in  current  was  observed  when  the  lithiation  potential  was  below 
3.42  V  and  delithiation  potential  was  above  3.435  V,  which  is  due 
to  the  phase  transformation.  The  CV  scans  starting  at  the  center 
of  potential  hysteresis  with  different  potential  amplitudes  demon¬ 
strated  solid  solution  behavior  inside  the  potential  hysteresis  and 
phase  transformation  characteristic  outside  of  the  potential  hys¬ 
teresis  range. 

EIS  has  been  widely  used  in  Li-ion  batteries  to  evaluate  the 
charge/discharge  kinetics.  A  reliable  EIS  can  be  obtained  only  if 
the  system  satisfies  the  criteria  of  linearity  and  time  invariance 
(LTI)  [23].  If  EIS  is  measured  at  the  open-circuit  after  charging  or 
discharging  to  a  given  state,  the  phase  transformation  occurs  only 
in  the  forward  potential  scan,  not  in  the  backward  potential  scan, 
resulting  in  a  SOD  or  SOC  shift  during  EIS  measurement.  Therefore, 
the  LiFeP04  is  not  a  linear  and  time  invariant  system  if  EIS  is  mea¬ 
sured  at  the  certain  charge  or  discharge  equilibrium  potentials.  The 
EIS  normally  measured  at  different  SOD  or  SOC  are  not  reliable.  A 
reliable  and  useful  EIS  can  only  be  obtained  by  shifting  the  testing 
potential  to  the  center  of  the  potential  hysteresis.  Fig.  10  shows 
the  EIS  plots  measured  at  50%  SOD,  50%  SOC  and  in  the  center  of 
the  potential  hysteresis  (marked  in  Fig.  4a)  with  different  poten¬ 
tial  amplitudes.  To  alleviate  OCP  shift  during  the  EIS  measurement, 
the  LiFeP04  sample  after  each  EIS  measurement  at  50%  SOD  (or 
50%  SOC)  was  fully  re-charged  (or  fully  re-discharged  for  SOC)  and 
then  discharged  to  50%  SOD  (or  charged  to  50%  SOC)  for  the  fol¬ 
lowing  EIS  test.  The  impedances  of  LiFeP04  in  Fig.  10  consist  of 
two  semicircles  and  a  sloped  line.  For  a  solid  solution  electrode,  the 
high  frequency  semicircle  in  EIS  is  attributed  to  the  SEI  film,  middle 
frequency  semicircle  to  charge  transfer,  and  low  frequency  line  to 
diffusion.  The  impedance  at  50%  SOC  is  much  smaller  than  that  at 
50%  SOD,  resulting  in  higher  rate  capability  in  charge  than  in  dis¬ 
charge.  EIS  in  both  50%  SOD  and  50%  SOC  are  less  sensitive  to  signal 
amplitude. 

Since  the  EIS  performed  at  charge  or  discharge  equilibrium 
potential  is  not  reliable,  EIS  should  be  performed  at  the  center  of 
potential  hysteresis.  When  the  EIS  test  is  conducted  at  the  cen- 
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ter  of  potential  hysteresis  using  amplitudes  less  than  half  of  the 
potential  hysteresis,  the  LiFeP04  can  be  treated  as  a  solid  solution 
since  no  phase  transformation  occurs  in  the  potential  hysteresis 
region.  When  the  amplitude  of  EIS  test  is  much  larger  than  half  of 
the  potential  hysteresis,  the  phase  transformation  can  happen.  In 
this  case,  if  the  kinetics  of  forward  phase  transformation  is  simi¬ 
lar  to  the  backward  phase  transformation,  a  reliable  EIS  can  still 
be  obtained.  The  EIS  measured  at  the  center  of  potential  hysteresis 
using  different  amplitudes  was  shown  in  Fig.  10.  The  impedances 
with  amplitudes  larger  than  half  of  the  potential  hysteresis  have 
a  larger  second  semicircle  than  those  with  amplitudes  less  than 
the  potential  hysteresis,  which  indicate  that  phase  transformations 
cause  additional  impedance.  The  phase  transformation  impedance 
appears  at  a  higher  frequency  than  Li-ion  diffusion  impedance 
and  is  mixed  with  the  impedance  of  the  charge  transfer  reaction. 
The  phase  transformation  impedance  of  LiFeP04  was  obtained  by 
subtracting  impedance  without  phase  transformation  from  phase 
transformation  impedance.  The  phase  transformation  impedance 
of  LiFeP04  is  around  20  £2.  The  phase  transformation  impedance 
is  50%  of  total  (phase  transformation  +  ion  diffusion)  impedance, 
which  is  similar  to  the  overpotential  ratio  (57%)  obtained  using  GITT 
in  Fig.  8. 

4.  Conclusions 

Equilibrium  potential  hysteresis  of  LiFeP04  is  a  thermodynamic 
property,  which  is  induced  by  strain  accommodation  energy.  The 
OCP  hysteresis  of  LiFeP04  in  the  phase  transformation  region  was 
determined  using  the  GITT  method,  and  the  corresponding  accom¬ 
modation  energies  were  calculated  according  to  the  potential 
hysteresis.  40  nm-LiFeP04  has  low  accommodation  energy  during 
lithiation  and  delithiation,  which  is  attributed  to  the  narrow  mis¬ 
cibility  gap  and  small  potential  hysteresis.  Phase  transformation 
only  occurs  at  a  potential  level  above  charge  OCP  or  below  dis¬ 
charge  OCP.  The  LiFeP04  inside  potential  hysteresis  region  behaves 
like  a  solid  solution.  Normal  EIS  obtained  at  different  SOD  or  SOC 
cannot  satisfy  the  criteria  of  linearity  and  time  invariance  (LTI).  A 


reliable  EIS  can  be  obtained  at  the  center  of  potential  hysteresis, 
and  phase  transformation  impedance  can  be  obtained  by  changing 
the  amplitude  of  the  EIS  signal  from  less  than  to  greater  than  half 
of  the  potential  hysteresis. 
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